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bJO' We experimentally study the shear flow of oil-in-water emulsion droplets in a thin sample chamber with a hopper shape. In 
^ .this thin chamber, the droplets are quasi-2D in shape. The sample is at an area fraction above jamming and forced to flow 
with a constant flux rate. Stresses applied to a droplet from its neighbors deform the droplet outline, and this deformation is 

quantified to provide an ad hoc measure of the stress. As the sample flows through the hopper we see large fluctuations of the 

^-H stress, similar in character to what has been seen in other flows of complex fluids. Periods of time with large decreases in stress 
lare correlated with bursts of elementary rearrangement events ("Tl events" where four droplets rearrange). More specifically, 
we see a local relationship between these observations: a Tl event decreases the inter-droplet forces up to 3 droplet diameters 
away from the event. This directly connects microscopic structural changes to macroscopic fluctuations, and confirms theoretical 
pictures of local rearrangements influencing nearby regions. These local rearrangements are an important means of reducing and 
redistributing stresses within a flowing material. 
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1 Introduction 
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[Newtonian fluids at low flow rates flow laminarly, with the 
'flow field independent of time. In contrast, complex mate- 

rials such as sand, foam, and emulsions behave like elastic 

solids at rest, and under sufficient applied stress, they flow. 
CO iWhile such materials have well defined time-averaged stresses 
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and velocity fields, their flows have strong stress fluctuations 
even at low flow rates This macroscopic observation is not 
surprising as microscopically the individual grains (or bub- 
bles or droplets) in a material must rearrange to allow for 
flow^, and their discrete size gives rise to these fluctuations. 
These rearrangements often occur in spatially heterogeneous 
■fashion—"—. Inter-particle forces are also spatially heteroge- 
]neous, with a small subset of particles bearing large forces, as 
'has been seen in experiments studying granular materials 
emulsion droplets^ii^, foams and simulations of friction- 
less particles ^ It is likely that small rearrangements of these 
[force-bearing particles influence stresses over a larger region. 

A variety of experiments have elaborated on this overall 
picture of local rearrangements leading to macroscopic stress 
fluctuations. In particular, flowing foams have been quite use- 
ful. Some experiments studied dry foams in 2D i^-is. ^jjggg 
are foams with bubble area fraction (|) > 0.95, where the bub- 
bles are deformed into polygonal shapes. Elastic stresses 
can be determined from the foam images, but this works 
only because of the polygonal shapes; this analysis is not 
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applicable to wetter foams (lower area fractions) where the 
bubbles have rounder shapes These dry foam experi- 

ments revealed many interesting behaviors, for example cor- 
relations between rearrangement events and gradients in the 
mean flow velocity and asymmetries between contracting 
flows and expanding flows Other experiments studied wet 
2D foams^ii^-i^. Due to imaging limitations, these experi- 
ments could not measure local stresses, only strain fields. Key 
observations from these experiments included observations of 
clusters of rearranging bubbles^ii^, some understanding of 
how the mean velocity field relates to the macroscopic rhe- 
ology''^'^""^^, and connections between instantaneous motions 
and time-averaged motions^. In general, these experiments 
focused on collections of rearrangement events and their re- 
lation to macroscopic stresses. Macroscopic stresses lead to 
average macroscopic strain profiles, which in turn cause col- 
lections of microscopic rearrangements that help the sample 
relax the macroscopic stress. The connection between indi- 
vidual local rearrangements and the relaxation of the stress 
has been difficult to see directly in experiments although con- 
ceptually it is clear such a link must exist. 

Motivated by these experiments, theoretical work attempted 
to connect these microscopic (particle-scale) forces and rear- 
rangements to the macroscopic stress fluctuations 'i^i^^i^^ and 
time-averaged velocity profiles Some theories-^Si^n^Z 

and simulations'^ make connections between individual mi- 
croscopic rearrangements, the local stress field, and the 
macroscopic flow. The general picture from these theories is 
that stress builds up locally, eventually exceeding a local yield 
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stress, leading to local rearrangements (perhaps in a cascade 
of several local events), which in turn relax the stress. These 
local events sum together to give macroscopic stress fluctua- 
tions. In particular, a specific prediction is that individual rear- 
rangement events relax the stress locally^^'^^ . This has not be 
directly tested in an experiment, due to the difficulty of mea- 
suring both bubble positions and stresses simultaneously. Un- 
derstanding the details of such localized rearrangements may 
help us understand phenomena such as shear-banding-?, and 
flow of complex fluids in general 

In this manuscript, we present experimental results of a 
"wet" oil-in-water emulsion where we can infer the stress 
on every oil droplet, and relate individual rearrangements to 
stress changes both locally and globally. Our experiment uses 
quasi-two-dimensional emulsion droplets. Small oil droplets 
are compressed between two parallel glass plates, deform- 
ing them into disks similar to experiments with photoelas- 
tic disks ^'^ or quasi-two-dimensional foams in Hele-Shaw 
j.gjjgjjji8j23,25j28 q^j. samples are jammed, and the stress 
each droplet feels is quantified by examining each droplet's 
deviation from a circular shape. In this experiment, we di- 
rectly and simultaneously observe macroscopic flow profiles, 
microscopic rearrangement events, and macroscopic and mi- 
croscopic stresses. In particular, we study microscopic "Tl 
events," where four droplets exchange neighbors, as shown in 
Fig- [na-c) ^'^'^i^'*i^^ '^^. During a Tl event, two neighboring 
droplets move apart and are no longer neighbors, while two 
other adjacent droplets move together and become neighbors: 
this is a change in topology. These rearrangements are induced 
by flowing the sample through a hopper [Fig. [Tld)]^i21^. We 
find that Tl events diminish the stresses felt by droplets over 
a distance up to ^ 3 diameters away from the event, demon- 
strating the existence of a flow cooperativity length scale, and 
showing how local rearrangements are responsible for macro- 
scopic stress fluctuations as predicted theoretically^^. 

2 Experimental methods 

Our droplets are silicon oil (poly-dimethylsiloxane, p=l g/mL, 
r|=350 mPas) in water, stabilized by Fairy™ soap of mass 
fraction 0.025, and are produced with the "co-flow" microflu- 
idic technique Our choice of emulsions (rather than foams) 
is for two main reasons: to avoid coarsening and to enable 
use of the microfluidic technique to control droplet sizes. Ad- 
ditionally, this is a new experimental system to complement 
prior experiments that studied granular media and foams. We 
place the droplets into quasi-two-dimensional hoppers shaped 
using a thin film of tape with thickness 0.10±0.02 mm. The 
tape is sandwiched between a glass coverslip and a glass 
slide. The walls are sufficiently clean that droplets do not 
have their contacts pinned to the wall, as can be seen by ob- 
serving the easy motion of droplets in a dilute sample when 




Fig. 1 (Color) (a)-(c) The three images show a typical Tl 
rearrangement of the four particles labeled by white circles. The 
field of view is 1.6x 1.2 mm", and the time interval between 
consecutive images is 0.66 s; images from Run 1. (d) Snapshot of 
sample. The field of view is 11.2x8.5 mm^. The flow direction is 
from left to right. The colored circles show the positions of Tl 
events taken from the time sequence of Fig.|2jb), where the earliest 
color (dark blue) corresponds to the Tl event at / = 38.7 s and the 
latest color (red) corresponds to the event at f = 47.3 s. Five of these 
Tl events occur in the time window t = (44.5 s,46.2 s) during 
which the droplets in that region move downstream only 0.5 mm; 
these five events are the leftmost events of the top group shown in 
the image. The background image, taken at ? = 45.2 s, has been 
inverted for better visualization. Each Tl event is positioned at the 
center of the four droplets comprising the Tl event based on their 
positions at f = 45.2 s. The data correspond to Run 1. (e) This 
image corresponds to (d), and is shaded based on the deformation D 
of each droplet. The darkest droplets correspond to D = 0.4 and 
white corresponds lo D — 0. 



the sample chamber is tilted. Moving droplets do experi- 
ence viscous forces from the walls (and each other), although 
these are generally weak enough that they do not deform the 
droplets directly In our experiment, these viscous forces 
are smaller than the inter-droplet forces—. This is because our 
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Table 1 Sample details of our six runs. The columns are flux rate A 
(mm^s^'), normalized flux rate A = A/7t{r)^ (s^'), area fraction (|), 
mean droplet radius (r) (mm), standard deviation a = \/JJj' 
normalized by (r), skewness s = (( ' J'^ )"^). hopper angle 
(degrees), and length scale ^/{r}- 



flow rates are small (leading to smaller viscous forces) and our 
area fractions are large (leading to larger inter-droplet forces). 
Nonetheless, larger droplets will certainly experience larger 
viscous forces from the walls, and likewise droplets moving 
faster will experience larger viscous forces. 

A constant flux rate is set by a syringe pump. This pump 
is attached to the far left edge of the flow chamber, which is 
^^50 mm to the left of the images of Fig.[nd,e). The chamber 
has parallel walls for 40 mm, and then begins to contract in 
the hopper region. It is this contracting portion imaged in the 
figure; note that the left region visible in the figure has paral- 
lel edges due to the limited field of view, but nonetheless the 
entire field of view is well within the contraction region of the 
chamber The mean velocity profile in this region is discussed 
below. We image the droplets with a 1.6x objective on an 
inverted microscope, using a 30 frame/s camera and a 0.33 x 
camera-mount for a large field of view. The flow rates are suf- 
ficiently slow compared to the camera rate that the droplets 
can be individually tracked using standard routines—. 

The details for our six experimental runs are given in Ta- 
ble[T] The samples have area fractions (|) > 0.90. At these high 
area fractions, all of our samples are in the jammed state^^i^. 
We find for our polydisperse samples that jamming occurs 
around (j)/ « 0.84: for (|) < (|)/, surface tension ensures the 
droplets are perfectly circular, whereas for our experiments 
conducted at (|) w 0.90, all of the droplets are slightly deformed 
even at rest. The standard deviation of the droplet radii is suf- 
ficiently large in all cases to frustrate long-range order Our 
results depend only on the mean radius (r) as will be described 
below, and do not otherwise change with different size distri- 
butions. 



3 Results and discussion 

3.1 Time-averaged velocity profiles 

Like the flow of other jammed material s '^i^i'^i^^i^^ ^ our time- 
averaged flow is smooth despite the complex motions of the 
individual droplets, with the velocity profile described by 



Run 


A 


A 





('•> 


a/{r) 


s 





X/{r) 


1 


2.93 


12.8 


0.90 


0.27 


0.21 


0.010 


25 


2.9 


2 


1.33 


18.8 


0.90 


0.15 


0.17 


-0.009 


25 


3.2 


3 


0.83 


5.5 


0.92 


0.22 


0.27 


-0.010 


27 


3.8 


4 


0.75 


14.1 


0.93 


0.13 


0.24 


-0.007 


27 


34 


5 


0.61 


9.9 


0.91 


0.14 


0.21 


-0.010 


26 


3.0 


6 


0.33 


6.2 


0.94 


0.13 


0.28 


-0.002 


27 


3.2 



(1) 



where y = is the centerline of the channel, a is the flow 
rate along the centerline, and P relates to the local strain rate. 
Droplets at the side walls slip along the wall with the velocity 
ofVx{x, > 0, where w{x) is the channel width that droplet 
centers can reach. The parameters a and p are proportional to 
the flux rate A as 



a = 



kaA 



(w(x)+2(r)) 



(w^(x) +6(r)w2(jc)) ■ 
1 .24 and ka 



(2) 



Equations[T]and|2l and parameters ka=l .24 and fep = 2.87 are 
all empirical observations. We find V,. ^ O.IV^ in all regions 
within the sample chamber The equations above ensurqj that 
the flux rate A is independent of position x. Our parabolic flow 
profiles are similar to those due to finite size effects seen in 
prior experiments, and likely reflect the finite size of our hop- 
per^". We emphasize that Eqn. [T] describes the time-average 
flow, and that the instantaneous velocity field can be and usu- 
ally is different. 

3.2 Global stress fluctuations vs. Tl frequency 

To quantify stress fluctuations in the flowing sample, we first 
examine the shapes of individual droplets to determine their 
stresses. Droplets are deformed away from perfect circles by 
forces from neighboring droplets^^. We quantify the defor- 
mation by determining the outline of each droplet, finding the 
radius r at each point on the outline (measured from the center 
of mass of the droplet), and then defining the droplet deforma- 
tion 

^' ' (3) 



the standard deviation of r normalized by that droplet's mean 
radius. D is determined for every droplet at every time, and 
in most of our results below we normalize by the mean value 
{D) within a region. Further details of calculating D are given 
in the Appendix. 

The physical interpretation of D relates to normal stress dif- 
ferences acting on a droplet due to neighboring droplets. More 

»■(■>-) 

* The total flux rate A (jc) is given byA{x) = / "^j",, V^{x,y)dy + 2V^{y = ■j)(r), 

where the second terni accounts for slip of droplets along the wall. Using 
Eqns. [T] and|2] A(x) = {ka — k^/12)A = A. That is, we only have one fit 
parameter, k^, as /:a = 1 +^p/12 to ensure constant (.v-independent) flux. 
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specifically, each droplet is acted on by compression forces 
from the surrounding droplets. If this compression is isotropic, 
then it contributes to the pressure on the droplet and does not 
change D. If the compression is not isotropic, there is a nor- 
mal stress difference acting across the droplet, and D is an 
ad hoc measure of this. Conceptually, one could consider the 
Cauchy stress tensor o acting on a droplet, and rotate the co- 
ordinate frame to find the principal normal stresses. In this 
coordinate frame the off-diagonal elements of the stress ten- 
sor are an ~ O21 =0. If there is a principal stress difference 
|c5ii — CT22I > 0, then this leads to an anisotropic deformation 
of the droplet into an ellipse-like shape. D is correlated with 
this principal stress difference. In a coordinate frame rotated 
by 45° from this one, the shear stress is maximal, and so D 
is likewise correlated with shear stresses; in fact, the maxi- 
mum shear stress is mathematically proportional to the princi- 
pal stress difference. Note that all nonzero deformation is due 
to forces between neighboring droplets. In experiments done 
at lower area fractions, below droplets do not touch one 
another and all droplets are completely circular (D ~Q within 
measurement error) at all flow rates. 

Figurellle) shows the spatial distribution of D at one instant 
in time. We find that the mean deformation obeys {D)yj (x) — 
Dq[1 +A/kvw{x)], with fitting parameters ky = 0.81 mm/s as 
a velocity scale and Do = 0.06 as the deformation for a non- 
flowing suspension at area fraction (|) « 0.9. In the absence of 
flow, {D)y — Dq; having {D)y > Do is because the nonzero flux 
rate results in droplets being deformed quicker than they can 
relax (limited by the viscosity). That is, our experiment is not 
in the quasi-static regime, otherwise (D) would not depend 
on the flux A. Nonetheless, the results we find below scale 
with the average deformation, and do not vary qualitatively 
from experiment to experiment. While our droplets are out of 
equilibrium, all relevant time scales (such as time scales for 
rearrangements) appear to be set by A for our data. 

The deformation rises near the walls and near the constric- 
tion, where the local shear rates are highest. In most of our 
analysis below, we focus on the left side of the sample cham- 
ber, where w{x) is large and thus (D)v,f is moderately indepen- 
dent of X. 

By averaging the deformation of all droplets within a large 
region as a function of time, we measure the global stress 
changes DG(f)i as shown in Fig. |2 a). These stress fluctuations 
are similar to those seen in granular hopper experiments ^i^^i^Q , 
In particular, the stress builds up and then can release during 
a short time interval, with the magnitude of the stress drop 
fairly significant (in many cases \ADdrop\ > 0.2(D)). While 
our data are insufficient to produce a clean power spectrum, 
we note that at high frequencies the spectrum is consistent 
with a power law P(co) ^ CO^^ This is true for fluctuations 
measured either in the left, middle, or right regions of our 
hoppers. This power-law decay is similar to various constant 
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Fig. 2 (Color online) (a) Temporal fluctuations of the mean global 
deformation Dq in the left shaded region of hopper as in the inset. 
This shaded region is ~ 8 mm in y-extent and so we use 
{D(j) = {D{w = 8mm)) = 0.087 to normalize the data. The data are 
from Run 1 (Table[Tll. (b) Subset of data from (a) indicated by the 
red dashed box, showing the fluctuations of the mean deformation 
within the left region of the sample chamber (c) The cumulative 
number of Tl events in that region over the same period of time. 
The vertical dashed lines in (b) and (c) indicate the time intervals 
when the mean deformation significantly drops, (d) The probability 
distribution of the frequency of Tl events fji divided by the 
normalized flux rate A=A/Tz{r)^, with the different curves 
corresponding to different regimes for deformation changes, as 
indicated in the legend. Open symbols correspond to deformation 
increases, and closed symbols to deformation drops. 



strain studies on the flow of granular materials i-^4,34- _36 yvijjcjj 
found P{(£)) ^ CO^' or ^ co^^, and is suggestive of similarity 
between our emulsion experiment (with only viscous friction) 
and these granular experiments (with static and dynamic fric- 
tion). 

These stress fluctuations are related to locaUzed rearrange- 
ment events, such as the Tl event shown in Fig.[Tfa-c). To il- 
lustrate this, in Fig.|2b) we plot a short segment of the global 
deformation Dcit) data from Fig. |2la). During this time pe- 
riod, we also count the number of Tl rearrangements that oc- 
cur in the same region, and show the cumulative number in 
Fig. |2c). Comparing these two graphs shows that Tl events 
happen more frequently during periods of large stress drops 
(36-38 s and 45-46 s). Furthermore, these events are spatially 
correlated, as shown in Fig. [Ud), which shows the locations 
of the Tl events occurring between 38 s and 48 s. Clearly one 



Tl event can trigger rearrangements of other nearby droplets^. 
(We count the Tl events using the algorithm given in the Ap- 
pendix.) 

We define a deformation drop as the magnitude of the de- 
crease of the deformation between a local maximum of D{t) 
and the subsequent local minimum (although first we smooth 
the data with a running average over a window of 0.33 s, to 
reduce maxima and minima which are only due to noise). 
Further evidence linking large deformation drops to groups 
of Tl events is found by calculating the frequency of events 
fri during a deformation drop. This frequency is defined as 
the number of Tl events occurring divided by the length of 
time over which the deformation decreases monotonically. We 
can similarly define fn during deformation increases. Prob- 
ability distributions of fn are shown in Fig. |2jd) for dif- 
ferent sized deformation changes. The solid triangles cor- 
respond to the largest deformation drops, and this distribu- 
tion shifts to higher frequencies. These distributions collapse 
across data sets when fn is scaled by the normalized flux rate 
A —A/n{r)^, and so Fig.|2d) shows data from all six runs. 
Figure |2d) shows that large stress relaxations are correlated 
with bursts of Tl events. Theory suggests that the rate of re- 
arrangements is a useful way to characterize the fluidity of 
jammed materials and our observations connect this fluid- 
ity to the large stress relaxations. 

3.3 Local stress relaxation around individual Tl event 

Given the correlation between Tl events and large stress drops 
(Fig.|2]i, we examine how a Tl event influences the deforma- 
tion of nearby (local) droplets Dl. We define Af = to be 
the instant of a Tl event [Fig.fljb)]. Distances from the event 
are measured by AR, where AR = is taken to be the cen- 
ter of the four droplets involved in the Tl event at time Af, 
that is, AR is measured in a co-moving reference frame. The 
mean local deformation Dl as a function of AJi is shown in 
Fig. [3ta) for several different Af 's both before and after the 
event. Before, the stress builds up at AR < 2{r), where (r) 
is the mean droplet radius; this stress buildup reflects strong 
deformation of the four droplets that will be involved in the 
Tl event. At the Tl event, the four droplets involved in the 
event dramatically decrease their deformation [solid circles in 
Fig. [3fa)], and quickly this stress release is propagated out- 
ward to distances AR = 6{r). To quantify this, we consider the 
quantity /(A/?) ==D(A/?,Afi) - D(A/?,Af2) using Afi = -0.5 s 
and Af2 = +0.5 s. We then compute 
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If /(A/?) — Aexp{~AR/X), then this calculation would yield 
the decay length X. For the six runs, we find X/{r) = 3.25 ± 



Fig. 3 (Color online) (a) Mean local deformation Dl around a Tl 
event, averaging over 186 Tl events during a 194 s duration movie 
(Run 1). These events all occur in the left side of the channel, see 
the inset to Fig.|2fa). The distance AR is defined in the frame of 
reference co-moving with the center of the four droplets undergoing 
the Tl rearrangement; for these data, this frame of reference moves 
roughly one mean radius per second. The Tl event occurs at 
Af = s. To reduce noise, the data are time-averaged over Af ± 0. 1 s. 
(b) Similar to (a), but as a function of Af for the given AR's as 
shown. Here the data are spatially averaged over a window of 
AR±0.2{r). For both panels, (Oc) = 0.087; see Fig. |2]for details. 



0.35 (mean and standard deviation). There is no obvious de- 
pendence on any of the experimental parameters, in particular 
the flux rate. The specific values of X for each run are listed in 
Table [T] Our directly observed length scale is comparable to 
those inferred from a three-dimensional experiment that stud- 
ied an emulsion with similar polydispersity^". Note that our 
data do not preclude the possibility of power-law decay, for 
which there would be no length scale, as has been theoreti- 
cally proposed"^. 

The temporal behavior around a Tl event is shown in 
Fig. |3jb), where the different curves correspond to different 
distances A/?. The black curve shows the changing deforma- 
tion for the droplets participating in the Tl event. Again, one 
sees a stress increase prior to the Tl event, a rapid drop at 
the event, then followed by a slower recovery after the event 
(likely limited by viscosity). The deformation drop seen in the 
dashed curve corresponding to A/? = 4(r) occurs about 0.5 s 
later in time, suggesting that the stress relaxation diffuses out- 
ward, as predicted by theory^*'. Together with Fig. [Tfd), the 
overall picture shown by our data is that Tl events cascade and 
release stress over a large region, leading to the fluctuations 
seen in Fig. EJa). One additional conclusion can be drawn 
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from Fig. |3] tlie mean local deformation is larger where a Tl 
event occurs. This is true even before the final increase right 
before the event; for example, for Af < — 3 s in Fig.[3jb), for 
the droplets closest to the future Tl event, {Di) « \.2{Dg)- 
The implication is that large stresses are often released at the 
point they are generated. 

Of course, Tl rearrangements have an inherent asymme- 
try between the diverging and converging droplet directions'^. 
To examine this more closely, we separately plot the behavior 
of Dl as a function of Af for these two directions in Fig. |4] 
Here it is evident that the diverging droplets show a much 
sharper decline in deformation after the Tl event. Droplets 
in those directions remain less stressed for quite some time 
after the event. In contrast, the droplets in the converging di- 
rection have a relatively rapid recovery of their deformation, 
approaching the mean value within a few seconds. (For this 
region of the sample chamber and this run, droplets move 
roughly one mean radius per second.) The combination of 
the data in these two directions yields Fig.[3jb), and it shows 
that the large stress release after Tl mainly comes from the 
diverging direction. This tendency is consistent among the 
six runs in Table [1] A similar spatial stress field was found 
in previous simulation'**, which did not examine the temporal 
behavior; Fig.|4]shows that this local stress field also changes 
with time. From Fig. |3]and|4] our observation confirms pre- 
dictions linking rearrangements to relaxation of stress over a 
larger region'^. Furthermore, Fig.S]highlights that such stress 
relaxation is not spatially isotropic, although this result is not 
surprising. Note that Figs. |3] and |4] show average tendencies; 
any given Tl event can raise the stress elsewhere, and this is 
likely one reason why clusters of Tl events occur [Fig. [lid)]. 
However, Fig. |3] makes it clear that on average, Tl events 
do indeed lower the stress of neighboring droplets and Fig. 2] 
demonstrates that this average tendency is true in both con- 
verging and diverging directions. The difference with prior 
simulations'^ is perhaps that the simulations considered the 
dry limit ((|) — > 1 .0) while we have a wetter system ((|) w 0.9). 

Our discussion has focused on Tl events, which are elemen- 
tary rearrangements within 2D samples. For 3D samples, Tl 
events would be replaced by more complex rearrangements, 
for example "shear-transformation zones" A more non- 
trivial difference is that our 2D experiment, every droplet feels 
a viscous drag force from the top and bottom walls, which 
would not be present for a flowing 3D sample Fortunately, 
as mentioned above the inter-droplet forces dominate over this 
viscous force. In particular, such viscous forces cannot cause 
rean^angements in 2D as they do not result in strain differ- 
ences between neighboring droplets. Thus, it seems plausible 
that our qualitative observations of local stress reductions re- 
sulting from 2D Tl rearrangements should still be relevant for 
3D flowing samples. 
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Fig. 4 (Color online) Similar to Fig|3jb), but we consider the local 
deformation £)l averaged along the directions of (a) the diverging 
droplet pair and (b) the converging pair. These directions are set by 
considering the instantaneous line joining the converging pair. 
Droplets within ±45° of that line (with the angle measured from the 
center of the four droplets) are along the converging direction, and 
the other droplets are along the diverging direction. In this way, 
space is divided into four equal quadrants, with the orientation of the 
quadrants changing with time as the Tl event evolves. For both 
panels, {Dq) = 0.087; see Fig.|2]for details. 

4 Conclusions 

We have used a quasi-two-dimensional emulsion to investigate 
the correlation between microscopic dynamics and macro- 
scopic stresses in a dense flow through a hopper Local re- 
arrangements ("Tl events") occur in bursts and are correlated 
with large stress releases. We observe a length scale for this 
correlation, where Tl events result in stress releases influenc- 
ing droplets as far as three droplet diameters away; this is 
the first direct observation of the "flow cooperativity length" 
predicted by theory'"'"' and suggested by previous experi- 
ments"^''^. Our results are for a system without static friction, 
and it is intriguing that the stress fluctuations we see are sim- 
ilar to those seen in granular experiments ^I'^i^^^^O fjjgsg simi- 
larities suggest that the connections we see between individual 
rearrangements, groups of these rearrangements, and macro- 
scopic stress fluctuations, may be common characteristics of 
complex fluids under shear for both frictionless and frictional 
systems. 

The overall implications of our results is that discrete mi- 
croscopic rearrangement events are an important means of 
reducing and redistributing stress within a flowing material. 
Clearly there is a strong connection between individual events, 
cascades of such events, and macroscopic stress fluctuations. 
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The localized nature of the stress reductions after a rearrange- 
ment event (Fig. |3]l show that stress fluctuations observed at 
the boundaries of a container are likely due primarily to re- 
arrangements near such boundaries. While flows of complex 
materials can often be usefully described by coarse-grained 
and time-averaged flow fields, it is clear that such averages 
hide interesting behavior: on shorter length and time scales the 
rearrangements are proceeded by significantly higher stresses, 
for example [Fig.|4fa)]. 
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Appendix 

We describe in more detail our procedures for determining the 
deformation D and the location of Tl events. 

5.1 Algorithm to compute deformation D 

First, a droplet is identified from its boundary in the micro- 
scope images, where the light is refracted at the interface be- 
tween oil and water. The raw images are grayscale images, 
and we determine a threshold intensity to convert the grayscale 
image to black and white: boundaries are black, and regions 
of pure oil or pure water are white. We then identify each 
connected white region in the image, and determine its area 
(total number of white pixels). Any concave regions are dis- 
carded as these are always voids. Additionally regions that are 
too small are discarded as either being the interstices between 
oil droplets, or else being oil droplets so small that no useful 
information can be gained from them. In practice, the area dis- 
tribution is bimodal and so it is straightforward to determine 
the cutoff. For the regions that are kept, we determine their 
centers of mass (giving all pixels comprising the region equal 
mass). To determine the perimeter, we consider rays drawn 
from the center of a white region outward at an angle, and 
identify the first black point encountered as a boundary point. 
For each droplet, we find 200 boundary points evenly spaced 
around the center (that is, spaced A9 = 27t/200 radians apart), 
and calculate their distance r(9) away from the center Using 
this r(9) data we can then calculate the mean radius (r) and 
deformation D for each droplet (see Eqn.|3)- 

The chief difficulty of very small droplets is our finite res- 
olution: the perimeter may only occupy a few pixels, and so 



r(9) is only coarsely sampled. One other problem arises when 
the droplet area fraction (|) approaches 1 . In this case, the voids 
between droplets become hard to see and this results in imper- 
fect identification of the droplet outlines at their comers. In 
fact, even to the eye the corners of the droplets appear to be 
distorted at area fractions above ^ 0.97. This is not a huge 
limitation; we expect that our algorithm for finding r(9) and 
thus D works quite well for area fractions up to ^ 0.97 and is 
less accurate for higher area fractions. 



5.2 Algorithm to identify Tl events 

One possibility for identifying nearest neighbor droplets is to 
use the Voronoi tessellation^"'^'. This partitions space into 
polygons, where each polygon consists of the points closer 
to the center of a given droplet than to any other droplet. 
Droplets whose Voronoi polygons share an edge are then con- 
sidered neighbors. This works well for monodisperse samples, 
but is less meaningful for our polydisperse samples. Accord- 
ingly, we use the Laguerre (radical Voronoi) tessellation in- 
stead'*""'*^. For this, when determine for a point which droplet 
center it is closest to, the Euclidean distance is weighted by the 
droplets' radii. This results in bigger droplets having larger 
polygons around them. Neighbors are still those droplets who 
share a polygonal face. This technique correctly identifies all 
touching droplets as neighbors. 

Knowing all neighbors at all times, the goal is to then iden- 
tify which neighbor pairs separate and which join together, 
and the times of these events. Unfortunately, the Laguerre tes- 
sellation (like the Voronoi tessellation) is sensitive to noise in 
the droplets' positions. For example, with a little bit of posi- 
tional noise, four identical droplets positioned at the corners 
of a square will flip between the two possibilities of diago- 
nally opposite droplet pairs being neighbors. To avoid over- 
sensitivity to noise, we require two droplets which are moving 
apart to not only cease being neighbors, but to also have their 
separation increase by at least 5% over a time interval of 1 s. 
Likewise, two droplets which are coming together must start 
being neighbors and also have their separation decrease by at 
least 5% over Af = 1 s. 

Having the information of the "moving apart" droplet pairs 
and the "coming together" droplet pairs, we search for events 
where we find two such droplet pairs. The final requirement 
is that the line segment joining the "moving apart" droplets 
should intersect the line segment joining the "coming to- 
gether" droplets. If this condition is met, then the intersection 
of those two line segments is taken as the position of the Tl 
event, and the time at which the topological change occurs is 
taken as the time of the Tl event. 
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